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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3283 


AERODYNAMIC FORCES, MOMENTS, AND STABILITY 
DERIVATIVES FOR SLENDER BODIES OF 
GENERAL CROSS SECTION 
By Alvin H. Sacks 

SUMMARY 


The problem of determining the total forces, moments, and stability 
derivatives for a slender body performing slow maneuvers in a compressible 
fluid is treated within the assumptions of slender-body theory. General 
expressions for the total forces (except drag) and moments are developed 
in terms of the geometry and motions of the airplane, and formulas for. 
the stability derivatives are derived in terms of the mapping functions 
of the cross sections. 

All components of the motion are treated simultaneously and second 
derivatives as well as first are obtained, with respect to both the 
motion components and their time rates of change. Coupling of the longi- 
tudinal and lateral motions is thus automatically included. A number of 
general relationships among the various stability derivatives are found 
which are independent of the configuration, so that, at most, only 35 
of a total of 325 first and second derivatives need be calculated 
directly. Calculations of stability derivatives are carried out for two 
triangular wings with camber and thickness, one with a blunt trailing 
edge, and for two wing-body combinations, one having a plane wing and 
vertical -fin. 

The influence on the stability derivatives of the squared terms in 
the pressure relation is demonstrated, and the apparent mass concept as 
applied to slender-body theory is discussed at some length in the light 
of the present analysis. It is shown that the stability derivatives can 
be calculated by apparent mass although the general expressions for the 
total forces and moments involve additional terms. 


INTRODUCTION 


Ever since R. T. Jones (ref. l) in 1946 demonstrated the use of 
Munk' s apparent mass concept of 1924 (ref. 2) for solving problems of 
slender wings in a compressible flow, an ever-increasing number of 
investigators have entered the field of analysis now commonly known as 
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slender-body theory. The stability derivatives of slender triangular 
wings were treated by Ribner (ref. 3 ) in 19^7 following the pattern of 
Jones, and in 19^8 Spreiter (ref. 4) extended the latter's result by 
means of conformal mapping to include certain wing-body combinations. 
Shortly thereafter, in 19%, Ward's general analysis for slender pointed 
bodies in steady supersonic flow (ref. 5 ) was published. 

After the appearance of Ward's analysis, a number of papers were 
written on various aspects of slender-body theory including extensions 
to subsonic flow and to "not-so-slender" bodies (e.g., refs. 6 and 7), 
and in 1952 Phythian (ref. 8) developed an analysis which included time 
variations in forward velocity and angles of incidence. Although many 
papers (e.g., refs. 9 and 10) have been devoted to the calculation of 
various stability derivatives for specific configurations, it is only in 
the past few months that a report by Miles (ref. ll) has given the com- 
plete counterpart of Ward's analysis for unsteady flow. 

The determination of stability derivatives has long been of concern 
to the engineer in connection with the dynamic behavior of airplanes, 
but the problem has assumed even greater proportions in the more recent 
slender configurations of missile design. The stability derivatives them 
selves correspond to the coefficients of a Taylor expansion representing 
a particular component of force (say lift) or moment as a function of 
the airplane motions. The coefficient of any particular motion (say q) 
in the expansion is equal to the partial derivative of the force or 
moment component with respect to that motion. Ordinarily, stability 
derivatives are defined as these partial derivatives evaluated with all 
of the independent variables except a set to zero, so that the usual 
stability derivatives depend upon the initial angle of attack as well as 
on the configuration. In the present paper, however, all derivatives 
are evaluated with all of the independent variables (a, 3, p, q, r, a, 3, 
p, q, r) set to zero. The advantages of this choice will become apparent 
in the course of the analysis. 

The present paper employs an approach believed to be novel in 
slender-body theory and is concerned with developing formulas for the 
forces and moments as well as the stability derivatives for general 
slender wing-body combinations. 1 The significance of the squared terms 
in the pressure relation for slender configurations precludes the 
possibility of considering the longitudinal and lateral motions inde- 
pendently, so all motions of the airplane are treated simultaneously. 


1 While the present analysis was being carried out, Bryson (ref. 12) 
published a paper treating essentially the same subject from a different 
viewpoint based on the tacit assumption that all the forces, moments, 
and stability derivatives can be obtained from the apparent mass analogue 
This assumption and some of Bryson's results are discussed in a later 
section. 
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The mathematical restrictions on the generality of the shapes that 
can be rigorously handled have been discussed in detail by Ward (ref. 5) 
and more recently for the unsteady case by Miles (ref. 11 ) . Such dis- 
cussion will not be repeated in this report. 
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LIST OF IMPORTANT SYMBOLS 


coefficient of — — term in series expansion of the mapping 
o n 


function t = f(cr) = cr + ~ 


n=o 


coefficient of — term in expansion of the complex potential 

p(£) 1 

value of A x at a = (3 = p = q = r = 0 

coefficient of Zn £ in expansion of F(C): B = — — 

2it dx 


distance from airplane nose to pivot -point 
complex potential cp + it 
length of airplane 

force in the z direction (approximately lift) 
rolling moment about the x axis 
reference length 

pitching moment about pivot point x = ci 
yawing moment about pivot point x = ci 
angular rolling velocity about the x axis 
pressure 

angular pitching velocity about the y axis 

fluid speed relative to axes fixed in the body 

component of q r normal to body contour in plane x = const, 
(positive into the fluid) 
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component of q r tangential to body contour in plane 
x = const, (positive counterclockwise looking upstream) 

V + iW 

angular yawing velocity about the z axis 

radius of transformed circle corresponding to airplane cross 
section 

cross-sectional area 
reference area 
time 

component of flight velocity along negative x axis 
component of flight velocity along positive y axis 
V Q - r(x-Ci) 

speed of a point fixed in the xyz system of axes 
components of q r in the x,y, z directions 
component of flight velocity along positive z axis 

W 0 - q(x-ci) 

force in the y direction (side force) 

Cartesian coordinates fixed in the body (x rearward, 
y to starboard, z upward) 

angle of attack (angle between arbitrarily chosen xy plane 
and flight direction) 

angle of sideslip (angle between xz plane and flight 
direction) 

angle between the positive y axis and the tangent to the 
body contour in a plane x = const. 

fluid mass density 

outward normal to the body contour in plane x = const, 
y + iz 
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complex coordinate of centroid of cross-sectional area 

(y c + iz c) 

complex coordinate in transformed circle plane 
velocity potential 

partial derivatives of 9 with respect to x,y, z, and t 

velocity potentials for unit velocity of the cross section in 
the y } z directions 

velocity potential associated with variations in shape and size 
of cross section with x 

stream function 

stream function along the contour of the cross section 


Special Notations 


contour integral taken once round the body cross section in the 
positive (counterclockwise) sense 


Force coefficients: Cy 


Moment coefficients: C m = 


(l/2) P U 0 2 S r 
M 


(l/2)pU 0 2 S r Z 


, etc. 

, etc , 


Stability derivatives: 



All derivatives are evaluated at 

* * 

a=(3 = p = q = r= a= p = 



3 2 C 



c^Cy 

CYpp = <5T 


= q = r 


= 0 


etc . 
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R real part 

I imaginary part 

9 

( ) derivative of ( ) with respect to time 

( ) complex conjugate of ( ) 

GENERAL ANALYSIS 


The problem to be treated here is the determination of the aero- 
dynamic forces and moments (except drag) and the stability derivatives 
for a smooth slender airplane or missile of arbitrary cross section per- 
forming slow maneuvers in a compressible fluid. The configuration will 
be limited in that the base (if any) of the fuselage and all wing trail- 
ing edges must lie in a plane essentially normal to the longitudinal body 
axis. 


Differential Equation and Pressure Relation 


The linearized differential equation for the velocity potential of 
unsteady motion of a compressible fluid is the well-known wave equation 



$ A A ‘ ®wi 


$ 




= 0 


( 1 ) 


where the system of axes Aq| is fixed relative to the undisturbed fluid, 
c G is the speed of sound in the undisturbed fluid, and t is time as 
measured in the Ap| system. Thus the velocity potential 0 is express- 
ible as $> = ®(A,p,|,t). 


In general, the pressure relation associated with the velocity 
potential $ is given by (ref. 13, p. 19 ) 

^ = - $ T - i qi 2 + J p x d + const. (2) 

where p 1 is pressure and q is the magnitude of the fluid velocity 
expressible as 2 

q i 2 = $ A 2 + V + (3) 


2 The subscripts on p and q are used to distinguish them from the 
angular velocities of rolling and pitching. 
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It will be convenient for 
the present problem to introduce 
a coordinate system xiyz which 
is fixed in the airplane. The 
axes chosen for this purpose are 
shown in the sketch and comprise 
a Cartesian coordinate system 
endowed with the translational 
velocities U Q , V Q , W 0 and the 
rotational velocities p, q, r 
of the airplane. (Wote that 
this does not constitute a com- 
pletely right-hand system.) 

The Xi axis passes through 
the airplane nose, and the ori- 
gin of the xiyz system is 
fixed at an arbitrary distance 
ci from the nose as shown in 
the sketch. 


( 



Since it is the purpose of this paper to study only instantaneous 
forces and moments (i.e., no time histories), it will be sufficient to 
choose an instant of time such that the positions of the moving xiyz 
system and the stationary system are just coincident. Thus, 

equations (l) and (2) will be expressed in the Xiyz system only for 
this instant, designated r = 0. For this purpose a new function 9 
is introduced such that 


c P(x 1 ,y,z,t) =®(A,h,|,t) (4) 

Wow, through the use of the transformations relating the moving and 
stationary coordinates (see e.g., ref. 13 , p. 12, and ref. l4, p. 79) 
one finds at t = t = 0 that 

$ = + ^9 8 x 1 + 89 8 y t 89 c)z 

T c)t 8 t dxi 8y St <3z 

= <Pt + u o<Pxi - (Vo + pz - rxi)9y - (W Q - py - qx x )9 z 

and 


= Vx x } ®n = % = ?z 

It can be seen from the sketch that the quantities (V 0 + pz - rxi) and 
(W 0 - py - qxi) are simply the velocity components in the y and z 
directions of a point fixed in the x x yz system. Note that in the 
corresponding x component (- U 0 + ry + qz) the products ry and qz 
are considered negligible compared with U Q . 
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With the assumptions of slender bodies, small angles, and slow 
maneuvers, the differential equation (l) reduces to Laplace's equation 
in planes xi = const, near the body; that is 

<Pyy + ^zz = 0 (5) 

It follows also that the density p must be treated as a constant 
in the pressure relation (2) which now becomes 

^ = - 3> t - | (f^ 2 + $| 2 ) + const. 

= - cp t - + (V 0 + pz - rx x ) 'Py + (W 0 - py - qxi)^ - 

i (<P y 2 + <P Z 2 ) + const. (6) 

Here, for the sake of convention, one can transfer the origin of the 
moving axes to the body nose by letting 


Xi = X - Ci 

so that pitching and yawing rotations are still made about an arbitrary 
pivot point x = ci* Thus, introducing the notation 

V = V 0 - rx x = V 0 - r(x - ex) 

W = W Q - qxx = W Q - q(x - Cx) 
equation (6) can be written 

^ - cp t _ u 0 9 x + (V + pz)T y + (¥ - py)9 z - i (9 y 2 + <P Z 2 ) + const. 

(7) 

This, then, is the pressure relation (referred to the moving body axes) 
upon which the calculations of the forces and moments will be based. It 
should be noted that a consistent application of the slenderness approx- 
imation requires the retention of the squared terms Ty 2 and Tg 2 . Thus 
slender-body theory is not a strictly linear theory although the differ- 
ential equation (5) is certainly linear. This means that solutions of 
equation (5) for T (and hence the velocities) can be obtained by super- 
position, but the pressures cannot. Likewise, the forces and moments 
cannot be calculated by superposition except for those special cases in 
which the contribution of the squared terms to the loading vanishes. 
Furthermore, when the airplane is performing combined maneuvers (e.g., 
simultaneous rolling and pitching) , the squared terms may contribute 
additional forces and moments. These in fact give rise to the second- 
order stability derivatives that will be included in the present analysis. 
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Total Forces and Moments 


The analysis to be presented here 
for calculating the total forces 
(except drag) and moments on a slender 
configuration will be the counterpart 
of a method originally due to H. Blasius 
(ref. 15 ) for obtaining the forces and 
moments on a two-dimensional body of 
arbitrary shape immersed in a steady 
incompressible stream. 3 This analysis, 
although not suited to the calculation 
of total drag, will nevertheless take 
proper account of the local forces 
associated with leading-edge suction. 
Consider a lamina of the slender air- 
plane, cut parallel to the yz plane, 
of thickness dx as shown in the sketch. 
One can write immediately the differ- 
ential lift and side force on an ele- 
mental area in terms of the local 
pressure p x on the body: 


z 



d 2 L = pi dy dx 

and 

d 2 Y = - px dz dx 


( 8 ) 


Wow, by introducing the complex variable £ = y + iz, one can express 
the differential complex force as 

d 2 Y - i d 2 L = - px dz dx - i pi dy dx = - i pi dx d[ ( 9 ) 


where £ is the complex conjugate of £ . In a similar fashion, the 
differential rolling moment about the x axis can be expressed as 

d 2 L' = - pi z dz dx - Pi y dy dx = - pi dx R(£d£) (10) 


where s denotes the real part. Further, the differential yawing and 
pitching moments about the pivot point x = c x are given by 

d 2 N - i d 2 M = - (d 2 Y - i d 2 L) (x - Ci) = i Pi(x - c x ) dx dl (ll) 

Integration of equations ( 9 ) > (10), and (ll) gives for the total forces 
and moments 


3 The method of Blasius has been extended to two-dimensional unsteady 
incompressible flows by L. M. Milne-Thomson (ref. l6) and recourse will 
be had to many of his techniques throughout the present analysis. 
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Y - iL = 


N - 1M = 


where the contour integrals are taken round the boundary of the airplane 
cross section in the positive (counterclockwise) sense. Before these 
integrations can be effected, the pressure p x must, of course, be 
expressed as a function of the complex variable £. Toward this end, 
it will be convenient to introduce two new definitions pertaining to 
velocities in the fixed and moving coordinate systems. First, the square 
of the speed of a point fixed in the xyz system can be written as 

Vi 2 = U Q 2 + (V + pz) 2 + (W - py) 2 (13) 

Second, it is noted that the square of the fluid speed relative to the 
xyz system is given by 

q r 2 = u r 2 + v r 2 + w r 2 = ( <P X + U 0 ) 2 + (<P y - V - pz) 2 + (cp z - W + py) 2 

(1^) 

so that, neglecting 9 X 2 in comparison with <Py 2 and C P Z 2 , one can write 
equation (j) in the form 

^ “ 9 t " \ !r 2 + \ V i 2 + const - (15) 

This expression will now be formed as a function of £ through the 
introduction of the complex variable R = V + iW and the complex 
potential F = 9 + i 1 l f . 



The speed Vx is immediately expressible as 


Vl 2 = U 0 2 + 


(V + pz) + i(W - py) 
= U 0 2 + (R - lp£) (R + ipD 


(V + pz) - i(W - py) 


(16) 


while the components v r and w r of the relative fluid velocity are 
related to F through the complex velocity by 
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Furthermore, at the body sur- 
face it will be seen from the 
sketch that 


+ iw r = q g e 


ie 


In 6 " 1 ^ - 0 ) 
= (q s - iq. n )e ie 


so that 


•v r - iw r = (q g + iq n )< 


-ie 


( 18 ) 



where q s and q n are the 
tangential and normal com- 
ponents of the transverse 

relative velocity and 0 is the angle defined .in the sketch, 
parison of equations ( 17 ) and (l8) gives 


Cora- 


ls + i( ln = ^ H " R - ( 19) 

from which 

l s 2 = (~ " R " ipl) e 2i0 - 2iq n ^ - R - ip 'f) e* 0 - q n 2 

It is now noted that (see sketch above) 

v r 2 + w r 2 = q g 2 + q n 2 = q/ - U Q 2 - 2U 0 9 X 

= (if " * " e2ie " 2i<ln (af " * “ ( 20 > 
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Finally, the boundary condition 
at the body surface requires that the 
normal component of the fluid velocity 
relative to the moving xyz system be 
equal to the normal component of the 
forward flight velocity. Let C be 
the contour of the airplane cross 
section in the plane x = const, and 
let C x represent the projection on 
that plane of the contour at x + dx 
( see sketch) . If V is the outward 
normal to C at any point, and dv 
is the distance between C and Ci 
measured along the normal, then the 
above boundary condition is, to the 
present order of accuracy. 


In 



( 21 ) 


Thus the pressure relation of equation ( 15) can now be written (for 
points on the airplane surface) in the desired form 



u o<P x “ £ 


d£ 


-\ 2 


- ip£ ) e 2i ® + iU c 




dy 

dx 


dF 


R 



i (R - ip£) (R + ipC) + const. 


( 22 ) 


Reduction of the integrals . - Before making use of equation (22) in 
writing the integrals for the forces and moments, it will be useful to 
notice from the sketch on page 11 that the differential distances on the 
body contour in planes x = const, are related by the angle 8 so that 
dy = ds cos Q and dz = ds sin 8 where ds is the differential arc 
length, positive counterclockwise . Hence, 

d£ = ds e ie ; d£ = ds e" ie ; d£ = d£ e " 2i0 ( 23 ) 


so that the first integral of equation ( 12 ) for the complex lateral force 
can be written, after expanding the squared term in the pressure relation, 
in the form 
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iL = ip J" dx J + ipU Q J dx jf <P x d£ + 

i 7 ; J dx jf (^| d£ + i | J' dx^ (R + ipt) 2 d£ 

ip / R dx ^ ~ d£ + PP y dxjTc^-d^ + 


pU Q dx jT ^ ^ ds - pU Q J' dx jT ~ (R + ip£)ds - 

i ~ dx (R - ip£) (R + ip£)d.£ (24) 


Note that the constant term in the pressure relation contributes nothing 
to the contour integrals of equation (12). 

The nine contour integrals of equation (24) can be divided into 
three types: (a) integrals that do not depend on the velocity potential; 

(b) integrals containing the real potential 9; and (c) integrals con- 
taining the complex potential F. The first type can be integrated at 
once and these will be dealt with first. The second type will be reduced 
to integrals of the third type by determining the stream function on the 
boundary, and the third type will then be handled by the method of 
residues . 


It is first noted that 


dv 

dx 


ds is simply the rate of change of 


— dv 

cross-sectional area S and that <p f — ds is the complex conjugate 

' dx 


of the rate of change (in the x direction) of the moment of cross- 
sectional area. Thus, one can write 


7^ (r + ip£)as = R ~ + ip ~ (s£ c ) 

dx dx dx 


(25) 


where £ c is the complex conjugate of the position of the centroid of 
area of the cross section. 

The other two integrals of equation (24) that do not depend on the 
velocity potential can be conveniently evaluated by the use of Stokes' 
theorem which can be stated in complex form as (see ref. l6, p. 130 ) 
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Ik 


or 



(26) 


where S is the area enclosed by the contour. Thus, using the first 
form of equation (26), one can write 


(R + ip£) d£ = 2i 



ipOaS = - UpS(R + ip£ c ) 


(27) 


where S is the airplane cross-sectional area. Similarly, from the 
second form of equation (26), 


j> (R - ip£) (R + ipC)<H = - 2pS(R + ip£ c ) (28) 

Thus, all the integrals of the first type discussed above have been 
evaluated. 

Before introducing the stream function for the evaluation of the 
first two integrals of equation (2^), it is well to note that the time 
differentiation can be taken outside the integral sign with no diffi- 
culty, but the x differentiation cannot since the contour of inte- 
gration is itself a function of x. It is shown in Appendix A that 

/^ aI ■ t / 941 - f c "’s Tte iy + 1 1 “v 

where C is the contour of integration round the airplane cross section 
and the surface of the airplane can be expressed either as 

z = Z x (x,y) 
or 

y = Yi(x,z) 

Corresponding to these expressions for the surface are the expressions 
for the slopes of the surface 


dZi _ dv/dx 
8x cos 0 

8 Yx _ dv/dx 

8x sin 0 


and 
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so that, recalling the relations 
of equations (23), one can write 


8 P 


W = ^ ^ + 


s (<p * + iq V )ds 

(29) 

Wow, in order to express the 

required integral cpd£ in terms 

of the complex potential through 
the relation 9 = F + it, the stream 
function on the boundary will be 
obtained from the boundary condition 
of equation (2l) . That is, the 

total outward normal of the fluid velocity in the plane x = const, 
given by the sum of q n and the normal velocity of a point on the 
boundary considered fixed in the xyz system. Hence, (see sketch), 



is 


= (V + pz)sin 0 - (W - py)cos 0 + U Q — (30) 

os dx 

and it is recalled that sin 6 = dz/ds and cos 0 = dy/ds. The sense 
of ds is indicated by the arrow along the contour. Thus, integrating 
along the contour, one finds that the stream function on the surface is 
given by 


to = Vz + - pz 2 - Wy + ■=■ py 2 + U c 


^ ds + G(x,t) 


where G(x,t) is an arbitrary function of x and t. Wow t s can also 
be expressed as a function of the complex variables £ and R by noting 
that 

I(R£) = Vz - w y = - ~ (H - rI) 

and 


y 2 + z 2 - U 
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whence 


V = - i (r5 - rQ + + U, 


Jn / — ds + G(xjt) 

0 ' dx 


(31) 


The integral J) 9d£ is now expressible in terms of the complex poten- 
tial since 9 = F + it and the integral f td£ can be evaluated from 


equation ( 31 ). That is, if one sets i = t D + where tx represents 
the integral term of equation ( 3 l) , then 


td£ = © t n d^ + p tidC = © d £ + [fx£L - <1 £dti (32) 


/ dV 

— ds or 
dx 

simply U n — and it is recalled that p t, — ds = — (s£ e ), so that 
dx J dx dx 

upon evaluating t Q dC by the second form of Stokes' theorem (eq. (26)) 


and noting that <j G(x,t)d£ = 0, one finds 


*d£ = -S(R + ip£ c ) + U 0 ^o — - U 0 — (s£ c ) 

dx dx 


(33) 


where £ 0 is the complex conjugate of the point at which the integra- 
tion was begun on the boundary. The final expression for Ji' cpd£ 
is therefore 


yW = Jmt + i jf +d£ =jT?d^ + iU D Co || - iS(R + ipC c ) - iU o“™( s U 

( 34 ) 

and the time derivative is 
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Of the remaining integrals in equation (24) there is one which still 
requires modification before the method of residues can be applied. 
Specifically, since F = F + 2it one can write 




+ 2idt) 


£dF + 2i 


Cdt 


and the second integral can be 
and noting that t taken once 
The resulting expression is 


integrated by parts using equation (33) 

dS 

round the contour has the value U D - — . 

dx 



£dF + 2iS(R + ipt ) + 2iU 0 A (SL) 
c dx 


(36) 


so that equation (24) for the complex force can finally' be written 
(after collection and cancellation of terms) in the form 


Y - iL = ip J' dx jT Fd£ + p J' S ^ dx + ipU G jT Fdt, - 

x=l 


ipU n A Fd£ - pU c 


J o 

x=o 
. p P l 


dx 


UqCq A . S(R + 2ip^ e ) - U Q A ( S £ c ) 

dx 


x=l 


a- 


- 2 j dx d£ - pU 0 — dx - ip / Rdx ^ dF 

pp f dx i £dF + ipp f S(R + ipC)dx 


(37) 


It should be pointed out that several terms have vanished by virtue of 
the fact that the x axis passes through the airplane nose. In par- 
ticular, note that S = £ c = 0 at x = 0. 

The complex potential .- Although equation (37) appears quite - 
unwieldy, all of the contour integrals are now in a form which admits of 
evaluation by the method of residues. For a body moving through still 
air, as in the present problem, all velocities vanish at infinity and the 
complex potential F can be expanded in a Laurent series of the form 


00 



A n U,t) 
t n 


D(x,t) 


F = B(x) Zn£ + 


n=i 


(38) 
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where B(x) is a source strength and the coefficients A n (x,t) give the 
intensities of all the higher order singularities representing the 
desired body shape and motions. This expansion applies only for large 
values of £ , but since there are no singularities outside the body, the 
contour integrals can be evaluated around a contour sufficiently large 
to insure the validity of the expansion. The arbitrary function D(x,t) 
is of no concern here since it can contribute nothing to the contour 
integrations. For the determination of drag, on the other hand, this 
function would be required. 

From equation (38), the derivative of F is 


so that 



dF = B V " 1 PAn 

d£ ~ t " L ^ n+i 

n=i 

00 / 00 

2B T n A n f V 

t 2 “ t L ^ l L 

n=x \ n=a 



Now it is seen that the residue of (dF/d£) 2 is zero since there is no 
l/£ term in the expansion. Therefore, 


Also 



2rti Res 


( 



= 0 


(39) 



|| d£ = 2iti Res(^ = 2itiB 


(4o) 


Similarly, if one writes the conjugate function 


F = B(x) ln£ + \ + D(x,t) 


n=i 


r 


(41) 


it follows that 


and 


£dF = 


jf Fd£ = 


i ^ dC 

d£ 


2rtiB£ n ~ 2niA± 


(42) 


2iti Res ( £ ^ ) = + 2jtiAx (43) 

/ 
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The coefficient B 


89 

integral ft — - ds 
8n 


(and, hence, B) can be evaluated by calculating the 
in two different ways. Thus 





dV 

dx 


ds = Ur 


dS 

dx 


and, by virtue of Gauss’ theorem. 


2 it 



89 

8r 


rid0 


r=ri 


2 it 


B 


o r i 


rxd0 = 2itB 


where ri is the radius of a circle enclosing the body cross section. 
Hence, it is seen that 


B = B = — — 
2 it dx 


(44) 


The final expression for the complex lateral force is obtained by 
using equations (39) to (44) in rewriting equation (37)* That is. 


Y - iL - 2itpU 0 Ai x _£ + pU G 


S(R + 2ip^ c ) + U G M (S^y 


8t 


dx 

<1 


x=l 


2itp / Ml. dx + p P S M dx + 2itipp / A x dx + 


8t 


ipp j S(R + ip£ c )dx 
'-'0 

For the case of steady straight flight 



(45) 


the complex 


force of equation (45) reduces to that given by Ward (ref. 5 ) • Although 
equation (45) applies to slender airplanes having cross sections of 
arbitrary shape, it is of interest that in a large number of practical 
cases, it is possible to choose the x axis so as to place the center 
of cross-sectional area always along the axis and thus to make £ c 
equal to zero. The simplest example would be an airplane having mirror 
symmetry of area about both the y and z axes. If the wings have no 
thickness, this places no restriction on the wings themselves with regard 
to number of wings, arrangement, dihedral, camber, etc. Equation (45) 
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as it stands gives the total side force and lift on a slender configura- 
tion having a nose at the forward end. However , it can also be used to 
give the contribution to these forces of a segment of the airplane lying 
between the planes x = Zi and x = l 2 by simply evaluating Y - iL for 
l = l x and Z = Z 2 and subtracting the results. 

If one considers the transformation of the arbitrary cross section 
in the £ plane to a circle in the a plane 

00 

5 = f(cr) = a + £ ( k6 ) 

n=o 

it can be shown (see Appendix B) that the dependence of the residue A x 
(and, hence, Ai) on the rolling velocity p is determined by the form 
of f(a) and therefore by the shape of the airplane cross section. It 
is found that Ai is independent of p if a Q vanishes and if f(u) 
contains either only odd negative powers of a (n odd) or else only even 
negative powers of o (n even) . It can be seen that if n is odd then 
f(-cr) = - f(cr) and, hence, the cross section has symmetry about two 
orthogonal axes. The statement can therefore be made that for airplanes 
having symmetry about both the y and z axes (no dihedral) one can 
determine the total complex force if he knows only the complex potential 
due to pure translation in the yz plane. 

An expression for the pitching and yawing moments can now be 
obtained from equation (12) by a procedure exactly parallel to that used 
in obtaining equation (^5)* Making use of the foregoing evaluations of 
the required contour integrals, one finds that the resulting expression 
is 


N 


iM = - 2itpU 0 / (x 


ci) ^ dx 
dx 




S(R + 2ip£ c ) 


+ 


_d_ 

dx 


u 0 — (sy 


dx - 2itp 


/ \ SAi 

(x - Cl ) r-i 

8t 


dx-p/° (x-ci)S — dx 

J n St 


2itipp 


l 

(x - ci)Aidx - ipp 



ci) (R + ip t c ) S dx 


( 47 ) 


The evaluation of the integrals for the rolling moment L' is 
somewhat different due to the additional t appearing in the expression 
of equation (12), so that the integrand of . some of the contour integrals 
appearing here will be nonanalytic in the variable of integration. This 
precludes a direct application of the method of residues. Such an 
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integral, which does, in fact, arise is 
instance it can be shown with the aid of 


I J) ££dF. However, in this 
equations (26) and (31) that 


/tig as (M» 

Evaluation of other required contour integrals by means of equations (26), 
(31 ), and (38) yields 



R 


(S) « • » 


(b9) 


£dF = - 2rtiA! 

and it can be shown in a manner parallel to that of Appendix A that 


B /^ a ^ a (s/ 9 Edf )- I /IS Us {50) 

so that the final expression for the rolling moment can be written in 
the form 


L' = pU 0 R d cp £d£ + pR J' dx ~ jf cp£d£ - 2jtpl J RA^dx 


x=l 

ppR f l S(R? C - 2RC c )dx - pU 0 I fffA ( s ^c)<^ 


(51) 


or, since R(U5) = \ a(U), 


L' = i pU 0 R f Fd(^) + - pR 
2 J 2 1 

x=l 


a dx A r Ed(CO - 2*pl P RAidx 

St J J 0 


l - d 


ppR /' l S(R£ c - 2R^ c )dx - pU 0 I / R “ (S^ c )dx 


(52) 
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In the case of configurations possessing mirror symmetry about both 
the y and z axes (no dihedral) , the general expression of equation (51) 
can be greatly simplified. If one writes for such a case 

9 = V9 2 + W9 3 + p9 4 + 9' 

where 92 and 93 represent pure translations along the y and z axes, 

94 is a pure rotation about the x axis, and 9' is associated with 
the variation of cross section with x, it can be shown from symmetry 
considerations that 


R J> 9 2 U£ = R 1 9 3 U5 = R 9' £d£ = 0 


R A 95d£ = pR A 9 a U.S 


and, since 9 4 is the potential due to pure rotation about the x axis, 
it follows that 


B Ui na£ r° 


Therefore, recalling that £ c = 0 for configurations having symmetry 
about the y and z axes, one finds that the general expression of equa- 
tion (51) reduces to 


= pU 0 pR <J 9 4 £cL£ - 2itpl / RAxdx 


and it has already been pointed out that for these cases Ax is inde- 
pendent of the rolling velocity p. Thus, for symmetrical configura- 
tions, the rolling moment has now been expressed as the sum of two inde- 
pendent parts, one due to pure rolling and one due to pure translation. 


STABILITY DERIVATIVES 


The specific maneuvers to be considered here will now be defined 
so that stability derivatives can be determined. Rotations are performed 
about the xyz (body) axes 4 at an attitude defined by the angles of 

4 Rolling about the wind axis can be treated as a special case by 
the proper choice of the arbitrary body axes. 
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attack and sideslip at t = 0, and 
accelerations are permitted in the 
fixed vertical and horizontal 
planes. In particular, if the 
angles of attack and sideslip at 
t = 0 are denoted by a Q and 3 0 , 
respectively, then the velocities 
V and ¥ at any time t are given 
by (see sketch) 

V = U 0 p 0 cos pt + 

U 0 a 0 sin pt - r(x - c x ) 


W = - U 0 a 0 cos pt + 

U 0 3 0 sin pt - q(x - c x ) 





From these expressions one finds that 


R = V - iW = (U 0 3 0 + iU 0 a 0 )e“ i P t - (r - iq) (x - d) 


and 


= - ip(U 0 3 0 + iU 0 a 0 )e-iP t + e-iP* (U 0 3 0 + iU 0 a 0 ) - 

(r - iq) (x - ci) - itp(U 0 3o + iU 0 a 0 )e“ i P t 

It will be noticed that in the above expressions, the velocity U Q (along 
the x axis) is considered constant. This means that pure pitching and 
yawing motions (q, r ^ 0) are performed at constant angles of attack and 
sideslip, so that for such maneuvers the airplane follows a curved flight 
path. Now, setting t = 0 in the above expressions, one finds 

R = U 0 3 + iU Q a - (r - iq) (x - c x ) 

and 

^ = U 0 3 + iU 0 cc - ip(U 0 p + iU Q a) - (r - iq) (x - c x ) (54) 
dt 

and these relations can be substituted directly into equations (45), (47), 
and (51) for the forces and moments. It will be noted that the subscript 
on a, and 3 has been dropped. This means that for the rolling case 
(p 0) a and 3 of equation (54) are not the time rates of change of the 
actual angles of attack and sideslip since a and 3 are measured in the 
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fixed vertical and horizontal planes. The alternative would have been 
to define <x and p # as > the actual rates of change of a and (3 (i.e., 
a = a Q - |3p and (3 = po + ap at t = 0) in which case the coupling of 
p with a and p would be obscured by the definition. This would lead, 
for example, to a nonzero value of aCm^ when the maneuver consists of 
pure rolling at an initial angle of sideslip. This seems undesirable. 

In seeking stability derivatives for the general problem under con- 
sideration, it will be advantageous to employ the transformation of the 
airplane cross section to the circle (see Appendix B) . In this way, it 
will be possible to carry out differentiations of the forces and moments 
explicitly with respect to the airplane motions and thus obtain the 
stability derivatives in terms of the transformation without specifying 
the complex potential. Thus, from Appendix B, 


Rax - Rr c 


2 , — 


IP a o r o + a i a o + 


a 2 a l a 3 a 2 


and therefore 


dAx _ - 8R 
8t 1 <3t 



ip( a o r o 2 + a i a o + • 


+ Al o 

( 55 ) 

• ) ( 56 ) 


where a x is the coefficient of the 1 /a term of the mapping function 
and r Q is the radius of the transformed circle. It_is recalled that, 
as shown in Appendix B, A lo is simply the value of A x at a = p = p = 

q = r = 0 and is therefore directly associated with the shape of the con- 
figuration and the choice of axes. 


The stability derivatives will be obtained by partial differentiation 
of the forces and moments with respect to each of the ten independent 
variables a, p, p, q, r, a, P, p, r and second derivatives will be 
included] that is, there will be derivatives of the types 



« « • 

where all derivatives are evaluated at a=p = p = q. = r = a= P = p = 
q = r = 0. The reason for this choice (which is not customary) will 
become more evident later, but it can be seen at once that all deriva- 
tives defined in this manner are constant for a given configuration and 
that there will be "cross derivatives" of the type which will 

show the mutual influence of the longitudinal and lateral motions. Thus 
the total rolling moment due to sideslip, for instance, will be express- 
ible as 
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C 7 * — C 7 + aC 7 + |BC 7 + pC 7 + • • . 

6 P 4 P ‘■a.p 4 PP 4 Pp 

For the sake of consistency, all coefficients and derivatives are based 
on the same reference area S r and reference length l r . 

If equations (U5) and (4-7) are rewritten with the aid of equations 
(54) to (56), the necessary partial differentiations can be effected with 
no difficulty and all the derivatives of the side force, lift, yawing 
moment, and pitching moment are obtained. The derivation of the rolling- 
moment derivatives is not quite so straightforward as the others since 
the expression of equation (51) for the rolling moment contains integrals 

of the type K, jT 9£d£ which cannot be handled directly by residues, 

as mentioned previously. However, when this integral is differentiated 
with respect to any motion except p or p, the resulting integral can be 
related to one of the integrals already evaluated by residues. In partic- 
ular, if we write once more 

9 = V9 2 + W9 3 + p9 4 + 9’ 

it follows from the boundary condition of equation (30) that 


S9 4 _ ^ dy + ^ dz 
8n ds ds 


or 


R(£&£) = 

Thus the integral appearing in 


y dy + z dz = 


89 4 

8n 


Cjp, for example, 


ds 

will be 


K /li ul = D °/' Pa lf d5 

But, by virtue of Green's Theorem (see ref. 13, p. 46), 



and again from equation (30) it can be seen that 
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8 c P g = dz 
8n ds 


so that, finally, one can write 


R 


^ US - u Q 


9 4 dz = - U 0 I 



This latter integral has already been evaluated with equations (34) and 
(te) to calculate C Yp . Similarly, the integrals appearing in the other 

rolling-moment derivatives can also be evaluated and expressed in terms 
of the mapping function of the cross section. The actual differentia- 
tions to obtain the rolling-moment derivatives are simple enough if one 
notes that the order of differentiation of the potential 9 is important 
in that the expression of equation ( 51 ) has been written for a specific 
instant of time (t = 0) . That is, 9 must be differentiated with respect 
to time first, then integrated to give the rolling moment, and finally 
differentiated with respect to the desired motion. Thus, since 9 is 
linear in the angles a and |3 as well as in the angular velocities 
p, q, and r, one observes that 

d 2 9 _ 8 2 9 

8a8t 8p8t 

and further that 

d g 9 89. 

8<x8t 8a 


8 2 9 


8 2 9 


8p8t 8a8p 


8 2 9 

8p8t 


89 . 

8p’ 


etc. 


The resulting expressions for the stability derivatives are given 
in table I, which is arranged so that all the side-force derivatives 
appear in the first column, all the lift derivatives in the second, and 
so on. It is found that a number of derivatives vanish identically, 
that is, regardless of the shape of the cross section. As a matter of 
fact, all but 84 of the possible 325 first and second derivatives vanish 
identically. For obvious reasons, the stability derivatives that are 
identically zero are not listed, but a definite pattern can be seen in 
table I which shows, for instance, that all second derivatives of C^, 
Cy, C m , and C n vanish except those involving p and that there are no 
second derivatives involving a, (3, p, q, or r. It should be noted that 
the order of differentiation is immaterial so that etc. 
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(if this were not the case, the total number of possible derivatives 
would be 550.) In the next section, a total of 49 relationships are 
found among the stability derivatives, so that, at most, only 35 deri- 
vatives need be calculated directly for any given configuration. It is 
important to note that the mapping function must be expanded in exactly 
the form of equation ( 46 ) before the proper coefficients can be obtained 
for use in the formulas of table I. In particular, the coefficient of 
the a term must be unity. 


Relationships Among the Stability Derivatives 


From the general results shown in table I, a number of interesting 
reciprocal relationships which are independent of the configuration are 
observed at once. For instance, the side force due to angle of attack 
is equal to minus the lift due to angle of sideslip (Cy a = - C;l ) . 

Similar equalities among the various derivatives are found throughout 
and are listed in table II. It should be noted at this point that many 
of the relationships of table II would be obscured by evaluating the deri- 
vatives at a ^ 0, as is customarily done. 


Beyond these simple equalities, there are some interesting relation- 
ships which can be brought out by an integration by parts. For example 
(see table i). 


c m a = " P l ( x " c i) ^-[2JtR(a! + r D 2 ) - S] dx 


S^»7- 


8x 


Syl 


r •'T 


(l - Cx)[2itR(a 1 + r Q 2 ) - S ] x=l - / [2itR(a 1 + r Q 2 ) - S]dx 


Ct + Ct . 

j_,q n a 


But from table I it is seen that 



Thus, one finds that 



( 


1 - ci 
l T 


) 
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This is a particularly interesting relationship in that it enables one 
to calculate from the static lift- and moment-curve slopes a quantity 
which would require dynamic tests in the wind tunnel for a direct exper- 
imental measurement. Another quantity of interest in this category can 
be obtained by integrating by parts the expression for C m given in 
table I. The resulting relation is " 



(58) 


which states that the damping in pitch is proportional to the lift-curve 
slope. This result, which is independent of configuration or choice of 
axes, was obtained previously by Bryson (ref. 12) whose analysis was 
implicitly restricted to bodies having a straight-line axis (i.e., no 
camber of the body or the wings). The apparent mass concept, which was 
the basis for the analysis of reference 12, will be discussed in a later 
section. Other relations obtained here in a manner similar to that for 
equations (57) and (58) are given in table II. 



It is of some interest to look into 
the damping in pitch of wing-body combi- 
nations on the basis of equation (58). 
Since the lift -curve slope is determined 
entirely by the trailing-edge configura- 
tion (see table i), it is evident that a 
wide variety of airplanes can be treated 
at once quite simply. If one considers, 
for instance, a configuration whose 
trailing-edge cross section consists of 
a circle with symmetrically placed 
straight lines, as shown in the sketch, 
it is known from the transformation (refs, 
k and 12) that the lift-curve slope is 
proportional to the quantity 



Thus, from equation (58), it follows that the ratio of the damping in 
pitch of the wing -body combination to that of the horizontal wing alone 
is also given by this quantity. It is important to note that any changes 
in shape ahead of the trailing edge (e.g., camber and thickness of the 
wings, variation in fuselage shape, etc.) are immaterial. Thus, for such 
configurations, the damping in pitch is plotted in figure 1 and it can 
be seen that (l) the body is always destabilizing, and that (2) this 
effect is a maximum for a body diameter to wing span ratio of l/JlT ■ 

The damping in pitch is made a maximum, on the other hand, by bringing 
the body to a point (or a line) at or ahead of the wing trailing edge. 
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The relationships given in table II are all independent of the con- • 
figuration, at least to the order of the present analysis. It can be 
seen from table I that for configurations involving any symmetries (for 
example if s£ c = 0), there will be additional relationships among the 
stability derivatives. Such a case will be considered in the section 
titled "Applications of the Theory." 


APPARENT MASS 


Perhaps the most striking feature of the results presented in table 
I is the frequent appearance of the quantities [2rtB,(ai - r Q 2 ) + S], 

[2rtH,(a;L + r Q 2 ) - S], and 2irl(ai). In fact, since both r Q 2 and S are 
real, the first two quantities are simply the real parts of 
[2rt(a! - r Q 2 ) + S] and [2it(ai + r G 2 ) - S], respectively, while the third 
can be written as the imaginary part of either of these bracketed quan- 
tities. It now becomes evident that a large number of the stability 
derivatives depend only on these two bracketed quantities which, in turn, 
depend only on the shape and size of the airplane cross section. As a 
matter of fact, it can be shown (ref. IT) that the quantities 
PR[2 rt(a x - r Q 2 ) + S] and pjRj^rt^a! + r Q 2 ) - S] are identical with the 
integrals defined in incompressible flow theory as the additional appar- 
ent mass of the cross section in the y and z directions, respectively. 


These are 


given by jT p^dz and jT p'Pgdy (refs. 13 and IT) where 


% 


and cpg are the velocity potentials for unit velocity of the cross sec- 
tion in the y and z directions. 


The mathematical basis for the use of the apparent mass to calculate 
the transverse force derivatives of slender bodies in steady flow (as done 
by Munk and Jones in refs. 1 and 2) was established with Ward's formula- 
tion (ref. 5) of the general expression 


Y 


iL = ipU G 


9d£ 


x=Z 


( 59 ) 


since a differentiation with respect to angle of attack, for example, 
yields 


d(Y - iL) 


3a 

X=l X=1 


Hence, taking the imaginary part of both sides, one finds 
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x=Z 


Uo 2 (m 3 3) x =l 


(60) 


where m 33 is the additional apparent mass of the cross section in the 
z direction for pure translation in the z direction. It thus becomes 
clear that the lift-curve slope is, in all cases, given by the apparent 
mass m 33 evaluated at the base of the body. If equation (59) were 
differentiated with respect to angle of sideslip rather than angle of 
attack, then the derivative Lp would be found to involve m 32 rather 
than m 33 , where m 32 is the apparent mass in the z direction for pure 
translation in the y direction. 


There remains now the question of the relation between the apparent 
masses and the total forces given by equation (59); that is, under what 
conditions can the total forces be calculated from the apparent masses? 
This can perhaps best be clarified by setting 


cp = v c cp 2 + W C <P 3 + 9 ' 


where the first two terms represent rigid-body translations of the cross 
section in the y and z directions and .9' represents variations in the 
shape and size of the cross section with x. For the steady case, the 
velocity components of the centroid of the cross section are given by 



and 



so that equation (59) can be written in the form 


iL = iU Q 2 


3 + Vm 32 - im 22 ) 

dx / 


a 


dz 


- )(m 33 - im 23 ) 


dx / 


+ 


J x=l 


ipU Q j cp'dC 
x-l 

Note that the quantities ! 3 + 


dye \ 
dx / 


and a 


dz c 

dx 


represent the angles 


that the line of centroids of the cross sections makes with the flight 
direction (i.e., the local angles of attack and sideslip). It can now 
be seen that the complex force of equation (59) is given by the apparent 
masses and the angles of attack and sideslip of the base cross section. 


provided that ft cp'd£ = 0. 
x=l 
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As an example, consider first a cambered body with no wings whose 
cross sections near the base are all circular. In this case 


q) 'x=Z = k ln 


J{j - y c ) 2 + (z - Z c ) 2 so that £ 9 ' d£ = 0 and the 


x=l 

apparent mass gives the side force and lift as well as their derivatives. 
On the other hand, consider a flat plate having deflected partial-span 
flaps. Here (when the plate is alined with the flight direction), the 
entire potential at the trailing edge is given by whose integral 


9'dC 


does not vanish. 


Hence, in this case, the lift-curve slope is 


x=Z 

given by the apparent mass but the lift is not. This is also true for a 
body of revolution having flat-plate wings at incidence to the body. 


For the unsteady case, as a result of a recently published report 
by Miles (ref. ll), one can show that the stability derivatives can also 
be obtained from apparent mass considerations. This is. most easily seen 
from the general expression of reference 11 


Y 


iL = 



9d£ 


+ 1 P 


dt 



<pa£ 


since, for example. 


(61) 



so that the imaginary part yields 


l 

m 33 dx (62) 


It should be noted here again that the total forces themselves are not, 
in general, given by the apparent masses. Miles also shows in reference 
11 that the rolling moment about the wind axis is given by 



L' = i pUoRjf Pd(SC) + | p ^ c M f ldx f ( 6 3) 

x=i 

Wow, by reasoning exactly parallel to that for the steady case, it can 
be concluded from equations (12), (6l), and (63) that in all cases all 
of the stability derivatives (except drag) for rigid slender bodies can 
be obtained from the apparent masses (or, more generally, the "inertia 
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coefficients"). Also obtainable from the inertia coefficients of the 
cross sections are (l) the total side force and lift, provided that 

cp'd£ = 0, (2) the total yawing and pitching moments, provided that 


x=l 


*>Z ~\ 

(x~Ci)dx — (J> cp'd£ =0, and (3) the total rolling moment, provided 

, 8x 


that K, 


cp'£d£ = 0. 


x=Z 


It is interesting to note that conditions (l) and 


(3) involve symmetries of the cross sections only near the base, but that 
condition (2) is more stringent and is satisfied, for example, by having 
bilateral symmetry of the cross sections over the entire length of the 
airplane . 


An alternative form for the integrals representing the apparent 
masses or inertia coefficients can be shown to be identical with the 
integrals representing the kinetic energy of the fluid associated with 
any desired unit velocity (linear or angular) of the cross section; for 

example, jf p9 3 dy = - jT p^Pg ds. However, it is essential to note 

that it is only for rigid-body motions of the cross section (as repre- 
sented, for example, by cp 3 ) that the two integrals are identical, since 

at the boundary. Thus, it can be 
even if T 1 is given proper 


^ CD y 

only for such motions does — = - — 

8n ds 


seen that jT p'P'dy £ - 


p cp ' ds 
Sn 


dimensions by dividing by a velocity. 


Inasmuch as relatively few inertia coefficients have previously been 
calculated, there seems to be little advantage (other than brevity) in 
expressing the stability derivatives in terms of these coefficients. It 
is felt that the formulas of table I involving the mapping function will, 
in general, be found more useful, although one should certainly make use 
of any of the coefficients already calculated. In this connection, the 
reader is referred to a recent paper of Kuerti, McFadden, and Shanks 
(ref. 17) in which the apparent masses of a number of interesting cross 
sections are listed for simple translation in the y and z directions. 


The apparent mass integral 


9 s dy was also calculated for a few 


shapes in connection with minimum drag problems in reference l8. How- 
ever, the integrations there were (for the rectangle) carried over both 
the exterior and the interior of the cross section since the configura- 
tion treated there was indeed a hollow rectangle made up of four thin 



NACA TN 3283 


33 


wings. ___ It is necessary then, in order to obtain the desired quantity 
[2rtH,(a 1 + r Q 2 ) - S], to subtract the cross-sectional area (as pointed 
out in ref. 17). It should be pointed out in this regard that all the 
results of the present paper are for solid bodies. Thus, if the body con- 
tains a jet, one must add to the calculated lift the negative rate of 
change of momentum (in the z direction) of the air passing through the 
jet. For a simple bent, circular, thin-walled pipe flying approximately 
along its axis, a jet exit velocity equal to the flight velocity just 
doubles the lift given for the solid circular cross section. 


The tabulated values of 


cPgdy in reference l8 have been adjusted 


for this internal flow and additional values have been calculated here to 
extend the range of the variable. The results are plotted in figure 2 to 
show the increase in lift-curve slope and damping in pitch obtained by 
use of a blunt trailing edge of rectangular cross section, of vertical 
end plates near the trailing edge, and of a biplane with sharp trailing 
edges. It will be recalled that both the lift-curve slope and the damp- 
ing in pitch depend only on the trailing-edge cross section, so the 
results of figure 2 are independent of wing thickness, camber, body shape, 
etc., ahead of the trailing edge. 


n fins 
n*3 


It might be mentioned 
here that if the apparent 
masses of a given cross sec- 
tion in two orthogonal 
directions are equal, then 
the apparent mass of the 
cross section is independent 
of its direction of trans- 
lation. This follows from 
the fact that the momentum 
vectors and the velocity 
vectors add in exactly the 
same fashion. In reference 
17 it was shown that the 
cross sections in the sketch 

possess this important property. It also follows, then, that ax = 0 so 
that many of the stability derivatives vanish for such configurations 
(see table i) . 




n corners 
regular polygon 


APPLICATIONS OF THE THEORY 


In this section, the results of the foregoing analysis will be 
applied to the calculation of the stability derivatives for several more 
or less special configurations. The first problems to be treated here 
will be concerned with the introduction of wing thickness and camber as 
parameters since the present analysis is applicable to unsymmetrical 
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configurations. Then, although the finned body of revolution has been 
treated by many authors (e.g., ref. 4, 9, 10, and 12), a plane wing-body 
combination will be considered in order to investigate the effect of the 
squared terms in the pressure relation. Finally, stability derivatives 
will be calculated for a wing -body -vertical-fin combination. 


Wings with Thickness and Camber 


The quantity "camber," as introduced here, may be complex, the 
imaginary part corresponding to the conventional camber in the vertical 
(xz) plane and the real part corresponding to a lateral camber or 
"wiggles" in the horizontal (xy) plane. Perhaps the simplest configura- 
tion of interest for the present problem can be made up of elliptic cross 
sections whose eccentricity and position in the lateral planes are arbi- 
trary functions of x. The required transformation for such a configura- 
tion is (see sketch) 



c 





a + 



(64) 


so that 


t a - b 

a o “ a i f j a 2 — 3-g — . . . — 0 

4 

qiK 

and the radius of the transformed circle is r D = . These quantities 

suffice for the calculation of many of the stability derivatives directly 
from table I. However, for the rolling derivatives, the complex potential 
is, in general, required. The complex potential in the transformed cir- 
cle plane can be derived from reference l6, page 239, and is given by 
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Fi(cr) = B Zn a 


H c (a + b)(b cos 7 + ia sin 7) 
2 a 


ip(a 2 - b 2 )(a + b)‘ 
16 a 2 


T n 


n=£ 


( 65 ) 


where H c and 7 define the lateral velocity of the centroid given by 


H c ei7 = R + U 0 


d|c 

dx 


iP^c 


(66) 


The logarithmic term of equation ( 65 ) permits a variation in size of the 
ellipse with x. The series term represents the portion of the potential 
required to meet the boundary condition of equation (2l) when dv/dx is 
arbitrary (warped body); this permits a variation of the eccentricity 
b/a with x. If now the potential of equation ( 65 ) is transformed to the 
t plane by means of equation (64) and the coefficient of the l/£ term 
is evaluated at R = p = 0, it is found that 



( 67 ) 

With this result, one can obtain all of the stability derivatives except 
Ci and Ci' directly from the formulas given in tables I and II, for 

any given configuration in this category. Two examples will now be 
considered. 


The first example will consist of a cambered elliptic cone; that is, 
an elliptic body with constant eccentricity whose span is a linear func- 
tion of x. The body axis will be chosen to pass through the center of 
the base and the camber line will be represented by a sine curve 
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sec. A- A 


z c = 8 sin rtx/c as shown in the 
sketch. This body can be thought 
of as a special case of a cam- 
bered triangular wing having a 
blunt trailing edge. The camber , 
of course, must be sufficiently 
small to insure the validity of 
the slenderness approximations. 

On the other hand, the thickness 
ratio as indicated by b/a is 
arbitrary. Thus the range of 
b/a from the flat plate 
(b/a = 0) to the circular cone 
(b/a = l) can be treated as one 
problem. It should be mentioned 
that the choice of body axes is 
arbitrary so that maneuvers about 
any other set of orthogonal axes 
fixed in the body could be 
handled equally well. A few of 
the interesting stability deri- 
vatives have been calculated for 
this configuration by the formu- 
las of tables I and II and the 


derivatives C n -n and C nm are 
ir CyP 

plotted in figure 3* It is seen that C n p is always negative for posi- 
tive camber, is linear in the camber, and increases with the thickness. 
On the other hand, Cn a p, which is negative in all cases, is independent 

of the camber and decreases to zero as the thickness ratio increases to 


one . 


As a second example of a configuration having elliptic cross sec- 
tions, the "wing-like" shape developed by Squire (ref. 19) is chosen. 

This shape has a variation of eccentricity of the ellipse such that all 
profile sections (except the midspan section) have a rounded leading edge 
and a pointed trailing edge. In addition to this particular thickness 
distribution, for the present problem the wing will be given a camber 
identical with that taken for the elliptic cone; that is, z c = 8 sin itx/c. 
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Thus, the wing to be treated 
here has the shape illustrated 
in the sketch. For this wing, 
the major and minor axes of 
the ellipse are given by 

a = ~ x (68) 

4 

and 

b = (c - x)x 
c 

where A is aspect ratio and 
t is the maximum thickness 
(thickness of the midspan sec- 
tion at x = c/2) . Some of 
the stability derivatives have 
been calculated for this wing 
by the formulas of tables I 
and II and C n p and C n(X p are 

plotted in figure 4. It is 
found that C n p displays, in 

general, the same variations 
with camber and thickness as 
did that for the elliptic cone . 



sec. A- A 


However, C n(X p displays a trend opposite 


to that for the elliptic cone. That is, C na p is seen to increase with 

thickness, so that the angle -of -attack contribution to the yawing moment 
due to rolling is apparently heavily influenced by whether the trailing 
edge is blunt or sharp. The derivative Cy^ was also calculated for 

both wings and was found to be independent of the thickness. In fact, 
for either wing, for the axes chosen. 



While this result appears to contradict the corresponding relation found 
in reference 12 (Cyp = Cjp), it simply highlights the fact that the analy- 
sis of reference 12 does not include camber although it could be extended 
to do so. Clearly, if 6=0 the two results are in agreement. The 
derivative Cj was also calculated for both cambered wings and was 

similarly found to be independent of the thickness. In fact, C i r was 
found to have a value equal to C n p for zero thickness; that is, 

C lr = - 0.4508 6/span. 


It is interesting to note that since the trailing-edge cross section 
of the Squire wing is a straight line, any stability derivatives that 
depend only on the mapping function of the trailing-edge section ( see 
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table i) will be the same as for a flat-plate wing (e.g., C T , C v , C T , 

L a *0 ijR 


C Lq> etc.). It is important also that for an elliptic cross section, the 


quantity [2jt(ax + r Q 2 ) - S] (apparent mass in the z direction) is inde- 
pendent of the eccentricity, depending only on the semimajor axis a, 
while the quantity [2jt(a! - r Q 2 ) + S] (apparent mass in the y direction) 
depends only on b. Thus it can be seen from table I that all of the 
purely longitudinal derivatives (C^, C Lq , C L& , C^, C ma , C mq , C m(i , C m ^) 

are the same as for a flat-plate wing of local semispan a. Similarly, 
all of the purely lateral derivatives (Cy , Cn^, etc.) are the same as 

for a flat-plate wing of semispan b if we replace (3 by minus a, Y 
by L, and N by M. One can see, therefore, that the ellipse is a very 
special cross section and tends to obscure some of the effects of thick- 
ness. For instance, it was seen in figure 2 that [2itH.(a 1 + r Q 2 ) - S] 
for a rectangle increases with the height of the rectangle and that a 
blunt trailing edge of rectangular cross section will therefore give an 
increase in lift-curve slope over a flat plate and a corresponding 
increase in the damping in pitch. 


For the evaluation of rolling moments it would appear from equations 
(5l) and (52) that some of the integrations might be quite difficult 
because of the nonanalytie character of the integrands. In fact the 
stability derivatives Cj p and C^ p contain the same nonanalytie inte- 
grands (see table I) . However, the integrations can sometimes be advan- 
tageously carried out in the transformed circle plane by the method of 
residues. For a configuration having an elliptical cross section at the 
trailing edge, for example, the calculation of the damping in roll C^ p 

becomes quite simple with this technique. Specifically, from table I, 


S r Z r 2 


R 


x=l 


^ a(«) = 

dp 


S r l r 2 


R 


x=l 


dF d(SS) 
dp da 


da 


(69) 


Now, since equations (64) and (65) for the required transformation and 
complex potential are already in the form of power series in a, and 
since on the circle boundary aa = r Q 2 , one can immediately write the 
integrand as a power series in a and therefore use the method of resi- 
dues. It is found for this case that the residue of this series is 


simply i (a 0 2 

8 



where the subscript refers to the trailing edge. 


Thus the damping in roll is given by 





TT 


R 


2rti 



flap 4 A _ bp 2 \ 

4s r l r 2 \ a Q 2 / 


(TO) 
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It is seen from equation (70) that for the Squire wing (b Q = 0) the damp- 
ing in roll is independent of the wing thickness and camber, remaining 
always at the flat-plate value as given by Ribner (ref. 3)* For the 
elliptic cone, on the other hand, C^p depends on the thickness and the 

effect of a blunt trailing edge becomes apparent. Values of Cjp for 

these cases are plotted in figure 5 . It should be noted in general that 
since 8 F/ 8 p is independent of any translational velocity in the cross 
plane, the damping in roll is not affected by camber. 

The device used above to permit the use of the method of residues 
can always be employed; that is, since 8 F/ 8 p, £, and £ are all express- 
ible as power series in a and a and since a and a are related by the 
radius of the transformed circle, the integrand becomes an analytic func- 
tion of the variable of integration in the transformed plane. However, 
if the transformation itself is an infinite series (as is the case for 
the finned body of revolution), then the residue and consequently the 
damping in roll will emerge as an infinite series involving combinations 
of all the coefficients of the transformation. This series is, in gen- 
eral, considerably more complicated than that entering into derivatives 
like Cy • 


Plane Wing-Body Combination 


It will be of some interest to consider here certain aspects of the 
plane wing-body problem in view of the fact that some stability deriva- 
tives had been calculated (see ref. 10 ) before it was generally realized 
that the squared terms in the pressure relation must be retained. A 
number of the simpler derivatives can be obtained quickly from table I 
if the mapping function of the cross section is known, and one of these 
will now be compared with the corresponding derivative obtained in ref- 
erence 10 without the squared terms in the pressure relation. The 
required mapping function is an infinite series obtained by making two 
successive Joukowsky transformations (see ref. 4), and it is found that 


a o — 



— £L^_ ~ 3- g — • • a — 0 

J 
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where s is the wing 
semi span and a is 
the body radius (see 
sketch) , Hence, many 
of the stability deri- 
vatives vanish by vir- 
tue of the fact that 
a x is real (see table 
I). It is noted that 
s, a, and £ c are here 
all arbitrary func- 
tions of x (within 
the slenderness 
approximation) and 

the body will have a straight axis (plane wing -body) only if £ c is a 
linear function of x. Also, if £ c is not a linear function of x 
then the evaluation of Ax (which is required for some derivatives) 

becomes a problem which amounts to determining Ax of the complex poten- 
tial for simple translation. This will be done shortly. 




For the pur- 
pose of illustrat- 
ing the influence 
of the squared 
terms in the pres- 
sure relation, the 
rolling moment due 
to sideslip will be 
calculated for the 
special case of a 
flat triangular 
wing mounted sym- 
metrically on a 
cone -cylinder as 
shown in the 
sketch. For this 


symmetry) there is no rolling moment provided by the portion of the body 
ahead of x = x Q and that the rolling moment due to sideslip at zero 
angle of attack is zero (Cjg = 0) . Therefore, the only pertinent deri- 
vative 5 to be calculated is which is given by (see table i) 



5 For this configuration all of the coupled (second) derivatives of 
the rolling moment vanish except C^g and Cjg^, and it is assumed here 

(for purposes of comparison with ref. 10) that q = 0. 
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(72) 


where S r and Z r have been taken as the gross wing area and the maximum 
wing span, respectively. The ratio of this to that for the wing 

alone (a Q /s 0 = 0) is plotted in figure 6 and is compared with the corre- 
sponding curve of reference 10. It can be seen that the error incurred 
by the omission of the squared terms in the pressure relation is in excess 
of 100 percent for ratios of body diameter to wing span greater than 0*5* 
Now since Cy a p = - C 1 according to table II, the side force due to 

rolling can also be compared with that obtained from a linear pressure 
relation by means of equation (72). This comparison is presented in 
figure 7 and it is seen that the difference is even more pronounced than 
that for the rolling moment due to sideslip. It should be mentioned that 
if the wing-body combination is cambered, the contribution of rolling 
moment by the nose will not, in general, vanish, nor will the rolling 
moment at zero angle of attack. 


It has been 
mentioned that for 
a general wing- 
body combination 
of the cross sec- 
tion discussed 
above (i.e., £ c 
not a linear func- 
tion of x) , the 
coefficient of the 
l/£ term in the 
complex potential 
must be determined 
if all of the stability derivatives are to be calculated. The coeffi- 
cient Ai for simple translation (p = 0) can perhaps most readily be 
obtained from the complex potential for a flat plate in a uniform stream 
by use of the transformation 




w = (5 - 5 C ) + -a — 

5 - Sc 

(see sketch). The complex potential in the w plane is given by (see 
e.g., ref. 16, p. l6l) 


Fi(u>) = - H c (w cos 7 - i sin 7 - d^) 


(73) 
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so that, by using the above transformation and adding a term H c £e“ i 7 to 
remove the free stream in the £ plane and a term B Zn ( £ - £ c ) to allow 
a variation of the radius a with x, one finds after expansion in series 
that the coefficient of l/£ is (for p = 0) 


Ai = 


H„ 


a 2 cos 7 


i (2a 2 - d 2 )sin 7 - b£, 


(74) 


where d = s + a 2 /s and it is recalled that Hce 1 -^ = R + U 0 (d£ c /dx) for 
pure translation. Now, by noting that for R = 0 


and 


cos 7 


Uo dy c \ 
H c \ dx ) 


sin 7 


Uo / dz c \ 
H c \ Hx ) 


one can write finally for this type of configuration 


A ln = 




- i (2a £ 
dx 2 


,2\ Hzq da 1 

d > dT + a S (y c 


+ iz c ) 


(75) 


With Ai q determined, all of the stability derivatives except and 


C^p can be obtained directly from tables I and II. 


The difficulties in 


determining and Cj* have been discussed in the preceding section. 


Wing -Body -Vertical-Fin Combination 



be exactly of 
reference 12 


In reference 12, the mapping 
function was developed for a body 
of revolution having four flat<- 
plate fins mounted 90° apart (see 
sketch) . Therefore, one can use 
the formulas of table I directly 
for such a configuration by first 
determining the proper coefficients 
in the expansion of the mapping 
function. It is important to note 
that for this purpose, the expan- 
sion of the mapping function must 
the form of equation (46). The mapping function given in 
is not of this form (as can be verified by carrying out the 
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expansion), but it can be modified to give the proper form. The resulting 
transformation is 



and it is found after a somewhat laborious expansion that the first six 
coefficients are 
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In general , 
the integrations 
indicated in 
table I would 
have to be per- 
formed numeri- 
cally for a given 
configuration of 
the type consid- 
ered above. How- 
ever ; for a con- 
ical configura- 
tion, great 
simplification is 
achieved and the 
integrations 
become trivial. 
Therefore, a two- 
parameter conical 

configuration is chosen to illustrate the effect of a vertical fin on a 
number of the stability derivatives; that is, the lower fin is removed by 
setting t 2 = a, and ti (or t) , a, and s are taken to be proportional 
to x (see sketch). Thus, the coefficients given above and the radius 
of the transformed circle become proportional to x or x 2 and most of 
the stability derivatives can be readily calculated. It is noted that 
for this configuration, with the axes chosen, A 1q = £ c = 0. Also, due to 

the conical property, C n p = = 0 if c x is chosen as 2/3 c. 

The stability derivatives obtainable from the coefficient a! (which 
in this case is real) and from the radius of the transformed circle r D 
have been calculated for a range of the parameters a/s and t/s and are 
plotted in figure 8. The purely longitudinal derivatives are of course 
unaffected by the vertical fin as seen in figure 8(b) (which is inci- 
dentally the same curve as given in figure 1 for a different purpose with 
regard to a more general configuration) . 


z 



-s 


Some interesting and important 
effects of the vertical fin can be 
seen in figures 8(a) and 8(c) which 
show a number of the lateral and 
"coupled” derivatives. First it 
should be noted that the two ends 
of the a/s scale correspond to 
the extreme configurations shown 
in the sketch. It is clear that 
the solid curves of figures 8(a) 
and 8(c) have no meaning for 

a/s > t/s since this would correspond to a vertical fin inside the body. 
Therefore, for values of t/s less than 1, the envelope curves (the 
dashed curves) corresponding to a/s = t/s have been plotted to fill in 


a/s=0 



a/s=LO 
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the rest of the range of a/s from 0 to 1. It is seen that even at 
a/s = 0 a vertical fin of small span is quite ineffective on both stabil- 
ity derivatives I and III (see figs. 8(a) and 8(c)). Hence, the ineffec- 
tiveness of small fins is not entirely caused by the body absorbing the 
fin, but is partly due to a blanketing effect of the horizontal wing. 

The significance of the effect of the vertical fin can perhaps best 
be appreciated by considering a relatively familiar derivative Cn a p 

(fig. 8(a)) which can be thought of as the angle -of -attack contribution 
to the yawing moment due to rolling. It appears that for fin heights 
such that t/s > 1.0 this derivative becomes positive over some portion 
of the a/s range (almost the entire range if t/s > 1.5). Thus, for 
this configuration, the angle -of -attack contribution to the yawing moment 
due to rolling can be changed from adverse (negative) to favorable 
(positive) by increasing the height of the vertical fin above about 1.5 
times the wing semispan, depending on the body radius. 

From the coefficients given in equation (77)* it is difficult to 
make any statement regarding the convergence of the series required for 
derivatives such as Cnp (see table i) . Therefore, one is at a loss to 

say how many terms of the expansion must be retained for satisfactory 
calculations. In order to get some idea of the convergence, the deri- 
vatives involving the series were calculated using 2, 3; 4, 5, and 6 
terms in the series and the results are presented in figure 9* It- appears 
that for some cases four terms would be sufficient. The results of fig- 
ure 9 indicate a strong blanketing effect of the wing on the vertical fin 
for t/s < 1. It can be seen that for a/t = 1 (plane wing -body) 

Cnp = C n p = . * . = 0. 

In reference 20, one of the configurations treated corresponds to 
the present conical wing-body -vertical-fin combination for a = 0 (no 
body) . As this furnishes an interesting check on the present calcula- 
tions, the appropriate values of Cyp, Cy^, and have been taken from 

that report 6 and are plotted on figures 8 and 9- It can be seen that the 
agreement is excellent, even for those derivatives calculated with only 
a few terms of the infinite series (fig. 9)* It should be noted in fig- 
ure 9 that wherever the best approximation curve (representing 6 terms 
in the series) cannot be seen, it is because the results were essentially 
identical with the previous approximation. 

It is interesting to note that for a cruciform wing-body combination 
(ti = ts = s), according to equations (76) and (77) and succeeding terms, 
one finds that a 0 = ai = a 2 = a 4 = a s = . . . = 0 so that many of the 
stability derivatives of table I vanish due to the symmetry of such a 
configuration if the axis of symmetry is chosen as the x axis; for 
example, here again Cnp = Cnp = • • . = 0. 

s The values of Cyp, Cyp, and C^p were taken from figures 11, 23, 
and 24 of reference 20 since There appear to be some typographical errors 
in equation (58) of that report. 
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CONCLUDING REMARKS 


A general analysis has been presented for determining the forces 
(except drag) and moments and the stability derivatives of a slowly 
maneuvering slender wing-body combination of arbitrary cross section. 

The results of the general analysis appear as (l) formulas for the forces 
and moments in terms of the airplane shape and motions and (2) formulas 
for the nonzero stability derivatives in terms of the mapping functions 
of the cross sections. 

Stability derivatives of the first and second order have been con- 
sidered so that the interdependence of the longitudinal and lateral 
motions is included. A number of relationships among the various stabil- 
ity derivatives were found which are independent of the shape of the air- 
plane, so that, at most, only 35 of a total of 325 first and second 
derivatives need be calculated directly. 

In order to bring out these relationships, the stability derivatives 
have been defined somewhat differently from the usual derivatives. For 
example, the usual derivative of rolling moment due to sideslip would, 
in the present analysis,- be given by 

C ^p = + “^ap + PC Zpp + pCz pp + q - Cj pq + rCz pr 

All derivatives as defined in this paper, then, are constant for a given 
airplane. Time rates of change of the angles and angular velocities have 
also been included, although these effects vanish for the particular 
derivative above. 

The use of the apparent mass concept for problems in slender-body 
theory has been discussed in the light of the present analysis and on the 
basis of previous treatments of slender-body problems by momentum methods. 
It is demonstrated that all of the stability derivatives can be calculated 
from the apparent masses (or inertia coefficients), but that the general 
expressions for the total forces and moments involve additional terms. 

From the results of the general analysis, some of the stability 
derivatives have been calculated for (l) two triangular wings having 
thickness and camber, (2) a plane wing-body combination, and (3) a wing- 
body -vertical-fin combination. These three cases have been used to show, 
respectively, (l) the effects of camber, thickness, and blunt trailing 
edge, (2) the influence of the squared terms in the pressure relation, 
and, (3) the effect of a vertical fin on the various stability 
derivatives. 

It was found that the effect of thickness on the angle -of -attack 
contribution to the yawing moment due to rolling was essentially opposite 
for blunt and sharp trailing edges, but the effect at zero angle of attack 
was similar in both cases. In both cases, the angle -of -attack 
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contribution was independent of the camber while the zero-angle-of -attack 
contribution was linear in the camber. 

The damping in pitch was shown to be proportional to the lift-curve 
slope and therefore to depend only on the base cross section. On this 
basis, the lift-curve slopes and damping in pitch were calculated for 
(l) a wing-body-vertical-fin combination, ( 2 ) a blunt trailing edge of 
rectangular cross section, (3) a sharp trailing edge with end plates, 
and (4) a biplane with sharp trailing edges. Sizable increases over the 
flat-plate values are shown in the last three eases. 

The derivatives usually called Cyp and were calculated to be 

more than 100 percent in error if the squared terms in the pressure 
relation are neglected in the case of a plane wing-body combination of 
body diameter to wing span ratio of 0.5 or greater. 

A number of stability derivatives were calculated for a conical wing- 
body-vertical -tail combination and the variations with body diameter to 
wing span ratio were plotted for various vertical tail spans. The influ- 
ence of the vertical fin was found to be markedly altered by the blanket- 
ing effect of the wing for small vertical tail spans. 


Ames Aeronautical Laboratory 

Rational Advisory Committee for Aeronautics 
Moffett Field, Calif., June 10, 1954 
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APPENDIX A 

DIFFERENTIATION OF A CONTOUR INTEGRAL 


z 



Consider a differentiation of the 


contour integral 



with respect 


to x, where the contour is a smooth 
closed path that depends on x, as shown 
in the sketch. The four points on the 
curve, designated a, b, c, and d, 
represent the maximum positive and 
negative values of y and z, so that 
the function 9 on the contour between 
a and c can be designated 9]_ ower (or 
simply that between c and a, 

9 U pper ( or ^u) > an( 3- those between b 
and d and between d and b, 9+ and 9-* 
respectively. Thus one can write 


_8_ 

8x 



9d£ 


f 9dy - i f 9dz = 

O Xup uX up OX 



4 - 


_8_ 

8x 



9 u dy - i ± 


d(x) 

b(x) 


9 + dz 



pb(x) 

'd(x) 


9_dz 


so that there are now four line integrals to be differentiated. This can 
be done directly by means of the formula 


8x 


is(x) 

■(x) 


f(x,y)dy = 


s (x) 8s „/ x 8r „/ x 

— dy + — f(s,x) - — f (r,x) 

r ( x ) ox ox ox 


which gives 
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A f Wt = r (X) 39lU<y ; Zl(X '. y) . 1 - dy + 1= cp (cx) - A 9j(a,x) * 

ox J 8x ox Ox 


r (x) a9u[x,y Zl (x,y )l 4y + |a . |= 

c ( x ) OX ox ox 


J f d(x) artfez,Ti(x,2)] Aa + aa cp + (d,x) -^9 + (t>,x)j- 

LJ b( x ) 8x 8x 8x j 


i r 9b(x) S9_ [x, z, Yi(x, z) l az + 3b „_ (bjx) . dd 9 _ (a>x) 
Ua(x) Bx Sx 


where z = Z x (x,y) and y = Yi(x,z) are alternative expressions for the 
contour C and must be single-valued between the prescribed limits of 
integration. 

Now, if the velocity potential <? is single-valued on the contour, 

then 

cp u = 9} at a(x) and c(x) 

and 

9 + = 9_ at b(x) and d(x) 

so that all four of the additional terms above cancel and one finds, after 
combining like integrals, that 


A/ W t = r ^tx.y.zxfe.yjl dy - 1 f MnnMnAI az 

8x J J c 8x J c 8x 


But the partial derivatives can be rewritten in the form 


^9[x,y,Z!(x,y)] [89 89 8 z x 


89 89 8 Z 1 

_ 8x 8z 8x _| r 


8x L8x 8z ox -l z =Z 1 (x,y) 


89[x,z,Yi(x,z)] [89 89 8 y x 


8x 8y 8x Jy=Y 1 (x,z) 
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Therefore 



, since 


89 

8x 


J z=Z x 


891 



the final result is 


9d£ 




89 8Z1 

8z 8x 


4y 





221Uaz 

8y 8x 



S9 8Z1 dy 
8z 8x 



Sl2hia z 

8y 8x 
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APPENDIX B 

THE RESIDUE A x OF THE COMPLEX POTENTIAL 


From equation ( 31 ) it follows that 


2i* a = r£ - R^ + ip££ + 2iU 0 J ds + 2iG(x,t) 

How, if one maps the body contour in the £ plane onto a circle of 
radius r Q , center the origin, in the a plane by the transformation 


5 = f ( a ) = a + y 5 


n=o 


(where a n is in general complex) then the "boundary function" of 
Milne -Thomson (ref. 16, p. 237) is obtained. Noting that on the circle 
boundary era = r Q 2 so that 


t - f (a) « — + 

a 


a n 
„ 2n 


a n 


n=o 


one can write the boundary function as 


2it s = R ( a + 




n=o 




ip a + 


n=o 


&n \ 

a n J 


r ° + V -5^ an ) + T (a) + 2iG(x,t) 


a 


n=o 


where T(a) is simply 2iU 0 


f 

— — ds expressed as a function of a. 
dx 


It has been demonstrated in reference 16 that this boundary function 
can be satisfied by setting the complex, potential equal to the part of 
2i 1 l f s containing only the negative powers of a. Thus one can set 


00 



n=i 



n=o m > n 


&m a n 

2n m-n 
r D a 
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where Tpj(cr) represents the portion of T containing negative powers of 
cr. The residue Ai of F can then be expressed as 


. . . ^ + A 1q 

where A 1q is the residue of T^-(a) and it is noted that A 1q can depend 

only on the shape of the cross section and not on any of the airplane 
motions; that is, A lo _must be a function of x alone and is simply the 
value of Ai at R = R = p = 0. Therefore A lo is zero for any con- 
figuration possessing an axis of symmetry if that axis is chosen as the 
x axis . 

It can be seen that A x is independent of the rolling velocity p 
provided that a Q vanishes and that either all the odd n or all the 
even n are absent in the expansion of f(cr). This leads to the con- 
siderations of symmetry given in the text. It should be noted that for 
symmetrical shapes a Q is the centroid of the cross section. 


Ai = Rai 


Rr f 


c 


+ ip + a i a o + 


a 2 a l a 3 a 2 
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TABLE I.- STABILITY DERIVATIVES 
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STABILITY DERIVATIVES - Concluded 
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TABLE II.- RELATIONSHIPS AMONG THE STABILITY DERIVATIVES 
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a 0 /s 0 


Figure!. - Damping -in- pitch for wing - body - vertical -fin combination. 
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(a) Yawing moment due to rolling at zero angle 

of attack , 


Figure 3 . - Effect of camber and thickness for blunt - 
trailing- edge triangular wing (elliptic cone). 
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(b) Angle -of- attack contribution to the yawing 
moment due to rolling or the rolling 
moment due to yawing . 


Figure 3~ Concluded . 



S r - wing area 
l r - span 



t/span 


(a) Yawing moment due to rolling at zero angle 
of attack . 


Figure 4.— Effect of camber and thickness for sharp- trailing - edge 
triangular wing ( Squire wing). 
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t/span 

( b ) Angle -of - attack contribution to the yawing 
moment due to rolling or the rolling 
moment due to yawing . 


Figure 4 .- C onduded . 




Figure 5 . - Effect of thickness on damping in roll 
for blunt - trailing - edge wing having elliptic 
cross section at the trailing edge . 
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Figure 6. - Influence of the squared terms in the 
pressure relation on calculations of the 
rolling moment due to sideslip . 




Figure 7 . - Influence of the squared terms in the 
pressure relation on calculations of the 
side force due to rolling . 
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Figure 8. - Stability derivatives for conical wing - body - 
vertical -fin configuration . 
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Figure 8- Continued . 



Stability Derivative 



a/s 

(c) Stability derivative ILL . 



Figure 8~ Concluded . 
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Figure 9 . - Calculation of stability derivatives for conical 
wing - body — vertical— fin combination (using different 
numbers of terms in infinite series of table I). 
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